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ABSTRACT 

Young OB associations with masses greater than 10^ Mq have been inferred to 
exist in the Galaxy but have largely evaded detection. Recently, a candidate OB 
association has been identified within the most luminous star forming complex 
in the Galaxy, the Dragonfish Nebula. We identify 18 young, massive stars with 
near-infrared spectroscopy from a sample of 50 members within the candidate 
OB association, including 15 0-type, and three Luminous Blue Variables or Wolf- 
Rayet stars. This number matches the expected yield of massive stars from the 
candidate association, confirming its existence and ability to power the parent 
star forming complex. These results demonstrate the existence of a 10^ Mq OB 
association, more powerful than any previously known in the Galaxy, comparable 
in mass only to Westerlund 1. Further, the results also validate the color selection 
method used to identify the association, adding credence to others discovered in 
the same way. 

Subject headings: infrared: stars — open clusters and associations: individual 
(Dragonfish Association) — stars: formation — stars: massive 



1. Introduction 



The most luminous and powerful of a galaxy's young stellar associations contain a signif- 
icant fraction of its massive stars, and dominate all coll ective forms of energetic feedback into 
the g alaxy's large-scale wind, gaseous disk, and halo (IMcKee fc Williams! 119971 : iDove et al. 
20001 ). While the massive cluster p opulations of some other galaxies, such as M83 are well 



characterized (jChandar et a. 



2010r). onlv the c losest couple kiloparsecs of the Milky Way 



have been probed in detail ( iLada fc Ladall2003l ). Across the Galactic plane, severe obscu- 



ration by dust and the high density of field stars have confounded searches for even the 
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most luminous associations, leaving our knowledge of them incomplete. About 50 young 
OB associations with masses greater th an 10"^ Mq have been inferred to exist in the Galaxy 
(IMcKee fc Williamsl 119971 : lLarsenll2009[ ): however, they have largely evaded detection. 



Our view of the G alactic population of H II regions is substantially more complete, 
thanks to radio surveys (ICaswell fc Hayned 119871 ) and subsequent follow-up studies. Of the 
identifiably coherent structures which emerge from such work, the most luminous is the so- 
called Dragonfish Nebula at (/,6) = (298,— 0.4) (IMurray fc Rahmarul2010l : iRahman fc Murray 
20101). Located on the far side of the Galaxy, 9.7 kpc through the Galactic disk, its morphol- 
ogy and radio velocities are indicative of an expanding bubble. Its H-ionizing luminosity 
of 10^^'* photons per second, if it originates from a coherent group of stars, requires the 
presence of the most luminous young stellar cluster or OB association in the Galaxy. The 
region's distance and location in the Galactic plane imply a high degree of extinction, roughly 
Ai^ - 1 or Ay > 10. 

Using a near-infrared (NIR) color selection method, a candid ate for the central pow ering 
source of the Nebula was identified: the Dragonfish Association (IRahman et al.ll201ll ). The 
association is composed of 406 ± 102 candidate O stars, consistent with the measured ionizing 
luminosity, making it one of the largest single OB associations or young clusters in the Galaxy 
(M ~ 10^ M0). Notably, the candidate association is not projected towards one of the bright 
regions of H II and 8 /xm emission in its vicinity, but within an evacuated shell surrounded 
by such emission. This suggests that the ionizing source has inflated a bubble of radius 69 
pc within the Galactic disk. The overdensity of sources being indirect evidence, we seek 
direct spectroscopic conflrmation that the overdensity consists of O stars, and that the fleld 
objects are K giants. 

In this Letter, we conflrm the existence of the Dragonflsh Association through NIR 
spectroscopy of massive stars. In addition, we identify two candidate Luminous Blue Vari- 
ables and a Wolf-Rayet within the association. This verifles the presence of a OB association 
which, on the basis of source counts, is sufficiently luminous to power its host star forming 
complex, giving it a total mass of 10^ Mq and making it the most luminous OB association in 
the Galaxy. The confirmation of the Dragonfish association validates the NIR color selection 
method as a way of identifying candidate clusters and associations in distant star forming 
complexes. 
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2. Observations & Data Reduction 

We conducted H- and K-band spectroscopy (1.52 - 2.4 /im) of 50 cluster sources pro- 
jected towards the candidate association, and a control sample of four non-cluster sources, 
between 2011 March 16 to 20. We use the SOFI spectrograph at the 3.6m ESO New Tech- 
nology Telescope with the red grism and 0'.'6 slit providing a spectral resolving power of R ~ 
1000. Both samples are limited to NIR colors consistent with extinguished O stars (or less 
extinguished K giants), and to K magnitudes brighter than 12. The control sample is offset 
from the candidate association by 0?5 and should represent the interloper sources within 
the cluster sample. The total integration time is 60 min for each cluster source and 2 min 
for each non-cluster source, using the nod-and-shuffle technique. This provides minimum 
signal-to- noise ratios of > 250 and > 50, respectively. The data was reduced and analyzed 
with Python, using the numpy and scipy packages. 



3. Analysis & Results 
3.1. Spectral Classification 



From the NIR color identification method, the contaminating non-member population 
confused wi th the candidate asso ciation is predominantly K giants with moderate extinction 
{Ak — 0.3; iRahman et al.l 120111 ). Due to the high extinction, the intervening interstellar 
medium will imprint absorption lines on the spectra as well, primarily from neutral metals 
(see §3.2p . The NIR spectra of K giants are similarly riddled with neutral metal lines, causing 
diffic ulty in identifying O stars on this basis, at least at low spectral resolution ( iRayner et al. 
2009h . 



To discriminate between extinguished O stars and K giants, we concentrate on two 
absorption lines which are characteristic of massive young stars but absent or very weak in 
K giants. These are the 2.166 //m Br 7 line and the 1.700 /um He I line, which is prominent in 
late O stars. A number of weaker features, such as the 1.693 /im He II line (present in early 
O stars), 1.681 /xm H I line, and the 2.188 fim He II line, are also useful for identification. 

In Table 1, we present the identified young massive stars: 15 O stars and three evolved 
massive stars. We present the spectra of a selection of these stars in Fig [TJ We classify 
the NIR spectra of O stars with the Brackett series of hydrogen absorption lines and He I 
and II absorption lines. These featur es are weak, necessitating signal-to-noise ratios above 
100 for detection (IHanson et al.l 120051 ). We determine a subtype range and luminosity class 
with the aid of the weaker absorption lines and NIR photometry using the NIR extinction 
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law from iNishiyama et al.l ( 120091 ) and model photometry from (iMartins fc Plezll2006l ). The 
magnitude differences between luminosity classes (from 0.2 to 1.0 mag in the K-band) of 
the same spectral subtype ranges make it improbable for small differences in the extinction 
determination to change the luminosity class identification. We note that massive stars with 
weaker features, such as early O stars lacking H e I absorption or wit h Brackett 7 absorption 
that has been filled in, may evade detection (IHanson et al.l 120051 ). All except the most 
luminous stars (the WN9 and 04-6 stars) show strong Br 7 absorption. The later stars all 
show He I absorption while the earlier stars show He II absorption. 

In additio n to the 15 O stars , we identify one Wolf-Rayet star based on H I, He I and 
He II emission (IMorris et al.lll996l ). Further, we find two candidate Luminous Blue Variables 
(LBVs), based on numerous He I emission lines, and [Fe II] emission in one case (Fig. |51 
Morris et al.lll996l ). These rare, exotic objects, transition stages of the most massive stars, 
are typically found in locations of recent star formation. Confirmation of the LBVs will 
require studies of their photometric variability. 



3.2. Interstellar Features 



The high signal-to-noise spectra of the massive stars presented in figure [T] contain a 
number of features that are similar to features in the contaminating K-giant population. 
Riddled throughout the spectra are neutral metal lines, primarily from Mg, Al, Fe, and Si. 
These elements are comr non in the interstellar med ium and have been observed in absorption 
within the local 100 pc ( iRedfield fc LinskyI 120041 ) . As such, these hues are expected to be 
imprinted on relatively featureless spectra when seen through a large column, as would be 
expected towards the association. 

In addition to the metal absorption lines, strong CO vibrational absorption bands are 
visible in the spectra of all identified massive stars at 2.29, 2.32 (z/ = 2 ^ 0) and 2.35 {ly 
= 3 — > 1) um. Th ese are often detected as photospheric absorption features of K-giants 
( iRayner et al.ll2009[ ). however, they may be also imprinted on nearly featureless spectra of 
younger, more massive stars by interstellar absorption. The ly = 2 — )■ bands have bee n 
detected towards NGC 2024 IRS 2 with line widths Av = 1.4 km s"^ felack fc Willnei]ll984h . 



The individual line s within each ba nd are found to be saturated with total band equivalent 
widths of ~0.5 A ( Lacy et al.lll994j ). Towards the Dragonfish region, there is CO emission 
detected between 10 km s~^ < visr < 30 km s~^ coinc ident with the Dragonfish associatio n 
and its surrounding photon-dominated region (PDR; iGrabelsky et al.lll988t lMurrayll201ll ). 
In the stacked cluster 0-star spectra (see §3.3p . the equivalent width of each of the CO bands 
is ~5 A, consistent with what is seen towards NGC 2024 when saturated over a larger range 
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of velocities. 



For the z/ = 3 — )■ 1 band to arise in the ISM, sufficient CO with excitation temperatures 
exceeding 1000 K must exist to populate the u = 1 vibrational state. The PDR traced 



by polycyclic aromatic hydr ocarbon emission in the reg ion ( iRahman fc Murrayll2010l ). is a 
likely site of such warm CO ( iHoUenbach fc Tielend 119971 ). We note that CO u = 2 ^ and 
z/ = 3 — > 1 bands have been observed in absorption towards massive stars in Cygnus 0B2 
classified by NIR and optical spect roscopy, at significantly lower dust extinctions than the 
Dragonfish's (jComeron et al.ll2008l ). Therefore, we conclude that CO absorption bands in 
the NIR spectra of the identified O stars are plausibly interstellar in origin, at least when 
the O stars are surrounded by an intensely illuminated PDR. However, this point deserves 
additional study. 



3.3. Stacked Spectra Comparison 



To investigate the weak NIR features in a more robust manner, we construct three 
groups of sources and examine the stacked spectra of each: the 15 cluster O stars, the 32 
cluster sources not immediately identified as massive stars, and the four non-cluster field 
sources from our control sample. We present the stacked spectra in figure [31 The primary 
spectral discriminators, Br 7 and 1.700 fim He I, are strong in the cluster O stars. Br 7 
is weaker but visible in the cluster contaminant population, and both lines are absent in 
the control non-cluster population. Further, weaker H I absorption features at 1.514, 1.534, 
1.544, 1.556, 1.570 and 1.681 /im, characteristic of massive stars rather than cool giants, 
are detected in the O star population and absent in the others. Conversely the Ca and Fe 
feature at 2.26 fim, a characteristic feature in the spectra of cool stars, is strong in the non- 
cluster population, weak in the cluster contaminants, and absent in the O stars. In addition, 
differences between the three populations are visible with lines at 1.522 /xm Mn, 1.529, 1.551, 
1.573, and 1.697 //m Fe, and 2.207 /im Na/Sc, all of which are present in the non-cluster 
and contaminant populations, and either weak or absent in the O star population. These 
spectral differences con&m that the identification of the association O stars is robust. 



3.4. Possible Runaway? 

One of the identified O stars. Source 14 in Table 1, appears to be highly blueshifted in 
the 1.693 /im He II, 1.700 fim He I, 1.736 /im H I and Br 7 absorption lines. The blueshift 
of the lines corresponds to a radial velocity of 300 ± 50 km s^^, as compared to the other 
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O stars that show no signficant Doppler shift. The velocity is much larger than can be 
explained by Galactic rotation. This indicates that the star is in a close binary system, or 
is a runaway O star either kicked out of the association or originating elsewhere. In our 
single-epoch spectrum, we find no indications of a companion star in the spectrum, such as 
splitting of spectral lines. If it is a runaway, it could have been produced by a binary-binary 
encounter or by the resu lt of the supernova explosion of the primary of a massive binary 
( iGvaramadze et al.ll2008i). Similar massi ve runaways have been identified originating from 
Cygnus 0B2 flComeron fc PasqualilboOTn . 



4. Discussion &: Summary 



The spectroscopic identification of massive, young stars confirms the existence of a 
young OB association within the Dragonfish Nebula, but its share of the nebula's ionization 
budget must be addre ssed statistically. Using the extracted candidate cluster properties from 
Rahman et al.l (120111 ). we infer a membership probability of 35.5 %, for an expected yield 
of 17.8 ± 3.4 members in the cluster sample. Indeed, 18 of our cluster sources are young, 
massive stars. This confirms that the Dragonfish Association is real, and not a chance sky 
overdensity of other stellar types (e.g. K-giants) coincident with the nebula or a spurious 
effect of intervening extinction. Further, this indicates that the NIR color selection method 
produces association membership probabilities consistent with the actual yield within the 
association. 

From the point sources identified usi ng the color selectio n method, assuming a min- 
imum candidate spectral type of 09. 5y (IMartins et al.l 120051 ) and using a standard stel- 
lar initial mass function (IKroupal 120011) . the inferred total stellar mass of the association 
IS ~ 10^ Mq fiRahman et al.l boilh . This is remarkably consistent with the estimated 
mass of the minimum stellar population required to power its central star forming com- 
plex (IRahman Sz Murrayll2010l ). This indicates that the Dragonfish Association dominates 
the ionization of its surrounding nebula. Further, the mass of this association is greater than 
other luminous associations in the Galaxy such as NGC 3603 (1.3 x 10^ M©), Trumpler 16 
in the Carina Nebula (1.8 x 10^ M©), Cygnus 0B2 (7.6 x 10^ M©), and the Arch es (7.7 x 10^ 



Clark et al. 



M^: IWeidner et all 120101 ). It is similar to that of Westerlund 1 (10^ M©; 
20051 ). however the Dragonfish is significantly more luminous compared to Westerlund 1 
based on its free-free lum inosity (Westerlund 1 is not visible in the WMAP free-free maps; 



Murrav fc Rahman|[201ol ) . This is consistent with the more evolved state of the massive stars 
in the latter. If the Dragonfish association is, in fact, the sole engine of the free-free emission 
in the complex, its output of 10^^'® H-ionizing photons per seconds is similar to R136, the 
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dominant powering source of th e most luminous H II region in the Local Gro up, 30 Doradus 
in the Large Magellanic Cloud (jParker fc Garmanyl Il993l : lEvans et al.l l201ll ) . 



The identification of O stars and evolved massive stars confirms the existence of the 
Dragonfish Association. The association is not embedded within a bright H II region, but 
rather a void of continuum and PAH emission surrounded by a shell, indicating that it is 
the central illuminating and infiating source of the Galaxy's most luminous star forming 
complex. As an extreme example of Milky Way star formation, the Dragonfish complex 
merits further study in a number of ways. Its large populatio n of O stars, an estimated 
13 with masses > 100 M0 (using the initial mass function from iKroupa] I2OOII ) . makes it an 
ideal laboratory for investigations of environmental infiuences at the upper limit of stellar 
masses. With high-precision astrometry, its membership can be more firmly established and 
its internal kinematics can be probed. X-ray, radio, and 7-ray follow-up is warranted to 
identify colliding wind binaries and young stellar objects. Searches for IR excess sources 
should reveal a population of new protostars triggered by its influence. Moreover, this 
influence should be sought in the dynamics of atomic and molecular gas in its environment. 

Addition al candidate associa t ions h ave been identified towards star forming complexes 
identified by iRahman &: Murrayl (120101 ) using this method and will be characterized in a 
forthcoming paper. The confirmation of the Dragonfish Association adds significant credence 
to the hypothesis that colour-selected clusters found toward the Galaxy's most luminous star 
forming complexes are in fact their driving stellar associations. 



This work is based on observations collected at the European Southern Observatory, 
Chile (ESQ 086.C-0502). We thank N. Murray, J. Graham, P.G. Martin and R. Breton 
for comments and discussion. D-S.M. and C.D.M acknowledge support from the Natural 
Science and Engineering Research Council of Canada. We acknowledge the staff of the ESO 
La Silla Observatory. This paper is dedicated to the memory of Crystal Marie Brasseur. 
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Fig. 1. — The spectra of 9 massive stars within the Dragonfish Association at the wavelength 
ranges of interest for spectral classification. The spectral types and identification numbers 
from Table 1 are indicated. The wavelengths of features used to classify the stars are indi- 
cated with dashed fines. The spectra have not been adjusted for Doppler shifts. The error 
bars are indicated in grey. 
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Fig. 2. — The spectra of the two candidate luminous blue variable candidates. The upper 
spectrum is star 9 and the lower is star 10 from Table 1. The error bars are indicated in 
grey. Characteristic spectral features are indicated. All spectral features are significantly 
stronger than those of the 0-stars. 
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Fig. 3. — The stacked cluster O star (top, black line), cluster contaminant (middle, lighter, 
line) and non-cluster (bottom, lightest line) spectra. The error of the stacked spectra is at 
the 0.01% level, and so is not indicated. Characteristic distinguishing spectral features are 
indicated for the O stars and K giants. The cluster contaminants likely include massive stars 
with features too weak to be identified on individual spectra. 
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Tabic 1. Massive Stars Identified in the Dragonfish Association 



# 


RA 


Dec 


1 (deg) 


b (deg) 


Spectral Type 


Kg (mag) 


Ak (mag) 


1 


12hllm30.2s 


-63°15'35'.'7 


298.458 


-0.736 


06-81 


10.5 


1.2 


2 


12hllm32.0s 


-63°15'52'.'7 


298.462 


-0.740 


08-9.51 


10.5 


1.1 


q 


1 OV. 1 1 m/1 1 Ac 

izniini4i.4s 


-Do iO oi.4 




n 7/1 s 
-U. ( 4o 


Uo-y.oiii 


ion 
iz.u 


1 n 
i.U 


4 


12hllm48.2s 


-63°14'33'.'6 


298.488 


-0.713 


06-8III 


11.1 


1.0 


5 


12hllm54.1s 


-63°17'04'.'0 


298.506 


-0.753 


WN9 


11.5 


1.2 


6 


12hllm58.4s 


-63°09'42'.'l 


298.495 


-0.630 


08-9.5III 


11.6 


0.9 


7 


12hl2m06.5s 


-63°11'36'/1 


298.515 


-0.659 


06-8III 


11.0 


1.1 


8 


12hl2mll.ls 


-63°12'2(y/8 


298.525 


-0.670 


04-6V 


11.6 


1.1 


9 


12hl2m20.3s 


-63°13'42'.'5 


298.546 


-0.690 


LBV? 


11.9 


1.1 


10 


12hl2m21.4s 


-63°16'38'.'4 


298.555 


-0.738 


LBV? 


9.0 


1.1 


11 


12hl2m22.1s 


-63°13'40'.'9 


298.549 


-0.689 


09-9.5III 


11.7 


1.0 


12 


12hl2m29.5s 


-63°10'51'.'5 


298.556 


-0.641 


05- 7V 


11.9 


1.1 


13 


12hl2m37.2s 


-63°14'52'/5 


298.580 


-0.705 


04-6V 


11.6 


1.1 


14 


12hl2m41.3s 


-63°09'27^3 


298.574 


-0.614 


<05V 


11.9 


1.1 


15 


12hl2m52.3s 


-63°17'03'.'9 


298.613 


-0.736 


07-9.5III 


11.2 


1.0 


16 


12hl2m52.6s 


-63°18'30'/8 


298.618 


-0.760 


06-8III 


11.4 


1.0 


17 


12hl2m54.4s 


-63°17'56'/5 


298.619 


-0.750 


06-81 


11.1 


1.2 


18 


12hl2m55.3s 


-63°06'4(y/6 


298.593 


-0.564 


09-9.5III 


11.3 


1.0 



